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Recent experiments using scanning tunneling microscopy show evidence for the formation of surface 
alloys of otherwise immiscible metals. Such is the case for Au deposited in Ni(110), where experiments 
by Pleth Nielsen et al. [Phys. Rev. Lett. 71, 754 (1993)] indicate that at low Au coverage ( <0.5 ML), Au 
atoms replace Ni atoms in the surface layer forming a surface alloy while the Ni atoms form islands on 
the surface. In this paper, we present results of a theoretical modeling of this phenomenon using the re- 
cently developed Bozzolo-Ferrante-Smith method for alloys. We provide results of an extensive analysis 
of the growth process that strongly support the conclusions drawn from the experiment: at very low 
coverages, there is a tendency for dimer formation on the overlayer, which later exchange positions with 
Ni atoms in the surface layer, thus accounting for the large number of substituted dimers. Ni island for- 
mation as well as other alternative short-range-order patterns are discussed. 


I. INTRODUCTION 

Recent experiments show clear evidence for surface al- 
loying of immiscible metals. 1 7 The atomic-resolution 
scanning-tunneling-microscopy (STM) experiments of Au 
atoms deposited on a Ni( 110) surface carried out by Pleth 
Nielsen and co-workers 1,2 show that in spite of the fact 
that there is a broad miscibility gap in the phase diagram 
of Ni-Au, Au atoms are seen to replace Ni atoms on the 
surface plane, forming a surface alloy. The experiments 
indicate the formation of a pattern of Au dimers in the 
surface plane, with the displaced Ni atoms located in 
chains along the closed-packed (cp) direction of Ni(110), 
following the fee stacking of the substrate. 1 The concept 
of a surface alloy for immiscible metals was corroborated 
by low-energy ion-scattering experiments performed by 
Boerma et a/. 3 In their work, they found evidence for 
the existence of Au atoms occupying near substitutional 
sites in the top Ni layer at low Au coverages. These ex- 
periments, added to the evidence advanced by the STM 
results and the effective-medium theory (EMT) calcula- 
tions verifying the conditions for surface alloy forma- 
tion, 1 provide enough motivation to perform larger scale 
simulations in order to investigate this growth mode. 

In this work, we address two fundamental issues relat- 
ed to this phenomenon: (1) the exploration of the 
equivalent to the ordered structures for the case of sur- 
face alloys, that is, the mixing patterns likely to form on 
the surface plane, and (2) the examination of the energet- 
ics of the surface alloying process of immiscible metals as 
a general phenomenon with the goal of establishing a cri- 
terion for the development of such type of growth. This 
would allow us to understand the features that distin- 
guish the Au-Ni case from other A-B mixtures that do 
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not form surface or bulk alloys under any conditions. 

For this study on alloy surface phenomena, we apply 
the Bozzolo-Ferrante-Smith (BFS) method for alloys, 8 a 
semiempirical technique which has had considerable suc- 
cess in previous applications. In Sec. II, we present a 
brief description of the BFS method, and in Sec. Ill, we 
apply BFS to the problem of surface alloying of Au-Ni. 
We summarize our results in Sec. IV. 

II. THE BFS METHOD 

Since its inception, the BFS method has been applied 
to a variety of problems, 8 starting with the basic analysis 
of bulk properties of solid solutions of fee and bcc binary 
alloys (heat of formation, lattice parameter, etc.) and 
more specific applications like the energetics of bimetallic 
tip-sample interactions in an atomic force microscope as 
well as Monte Carlo simulations of the temperature 
dependence of surface segregation profiles in Cu-Ni al- 
loys. Other applications include surface structure of me- 
tallic alloys and a diagrammatic analysis of ordered alloy 
clusters for the determination of the ground-state struc- 
ture of a given binary alloy. An additional advantage of 
BFS is that it allows for deriving simple, approximate ex- 
pressions which describe the trends in segregation as well 
as elucidating the driving mechanisms for these phenom- 
ena. Also, as a consequence of the ideas underlying the 
foundation of BFS, simple expressions for predicting the 
composition dependence of bulk alloy properties based 
solely on pure component properties have been recently 
derived, providing an alternative to the commonly used 
Vegard's law. 

In what follows, we present a brief review of the 
method. Due to its way of partitioning the energy in 
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different contributions, this presentation should be com- 
plemented with a review of previous applications, in or- 
der to familiarize the reader with the main concepts dis- 
cussed below. 

The BFS method is based on the idea that the energy of 
formation of an arbitrary alloy structure is the superposi- 
tion of individual contributions e, of nonequivalent atoms 
in the alloy, 8,9 

e,=E?+g,(ef-ef 0 ) , (1) 


feet crystal, as shown below. Within the framework of 
ECT, U af is interpreted as the lattice parameter of an 
ideal, perfect crystal (i.e., the equivalent crystal), where 
the energy per atom is the same as the energy of atom i in 
the actual, defect crystal. 

In general, the ECT equation for computing the st ain 
energy reads 


NR p x e 


ctR 


1 + MR \e 


[a-Hl/X)]*. 


= 2/; 


~[a + S(rj)]rj 


( 8 ) 


so that the total energy of formation is 


A// = 2e, . 


( 2 ) 


For each atom, we break up the energy into two parts; 
a strain energy e 5 and a chemical contribution, linked by 
a coupling factor g: 

e,- = e- ,rain +g J ef' em , (3) 

c 

where / denotes the atomic species of a given atom (e 0 is 
a reference energy to be defined later). 

The strain energy ef accounts for the actual geometri- 
cal distribution of the atoms surrounding atom i, comput- 
ed as if all its neighbors were of the same atomic species 
as atom i. ef is then evaluated with any available tech- 
nique. 

The coupling term g, is related to the strain energy in 
the sense that it contains information on the structural 
defect included in ef. In order to establish this connec- 
tion, based on the assumption that the universal binding- 
energy relationship of Rose et al 10 contains all the 
relevant information concerning a given single- 
component system, we write 

zf=EiF*(af*) , (4) 


where 


F*(a*)—\—(l+a*)e~ a 

and where af* given by 

(n, s -nj) 


*, 5 *=<r 


/: 


(5) 


( 6 ) 


is a scaled lattice-parameter related to af 9 a quantity that 
contains the structural information of the defect crystal. 
a f /,, and E l c are the equilibrium lattice parameter, scal- 
ing length, and cohesive energy of a pure crystal of 
species / and ^ 3 = (3/167r) for fee metals. 

Once ef is evaluated by any theoretical means, a?* can 
be easily obtained from Eq. (4) with which the coupling 
term g, becomes 


As in previous efforts, 8 we choose the equivalent crys- 
tal theory (ECT) (Ref. 1 1 ) to perform strain-energy calcu- 
lations, the choice being guided by the simplicity and reli- 
ability of this technique. Using ECT for computing ef in- 
troduces the added advantage that af (and thus a?* ) is 
directly obtained by solving the ECT equation for the de- 


(see Ref. 1 1 for details), where the quantities p , a, X and 
the screening function S are defined in Ref. 11. The <um 
on the right-hand side of Eq. (8) runs over all neighbors 
of atom i at a distance r ; . Equation (8) is then solved for 
the lattice parameter of the equivalent crystal af. R , and 
R 2 are the corresponding nearest- and next-nearest- 
neighbor distances in the equivalent crystal. The strain 
energy is then computed with Eq. (4). For the particular 
case where all the neighboring atoms are located at lat- 
tice sites, rj = r { and S(rj) = 0 for nearest neighbars, 
r j =r 2 and S(r 2 )= l/X for next-nearest neighbors and, if 
n is the actual number of nearest neighbors and m is the 
corresponding number of next-nearest neighbors, then 
Eq. (8) is simply 

NR-,e + 

= . (91 


Rigorously, the computation of the strain energy in- 
cludes four terms (see Ref. 11). In this work, we neglect 
the three- and four-body terms dealing with the bond an- 
gle and face-diagonal anisotropies and retain only the 
two-body term that accounts for bond-length anisotro- 
pies, 11 which we expect to be relevant for atoms in the 
top (surface) layers. The higher-order terms would be 
proportional to the small local fluctuations of the atomic 
positions around the equilibrium lattice sites. We expect 
that the leading term, Eq. (4), will adequately account for 
these small distortions. 

The chemical contribution ef is obtained by an ECT- 
like calculation. As opposed to the strain-energy term, 
the surrounding atoms retain their chemical identity , but 
are forced to be in equilibrium lattice sites. If N ik (M ik ) 
denotes the number of nearest (next)-neighbors of species 
k of the atom in question (of species /), then the ECT 
equation 11 to be solved for the equivalent lattice parame- 
ter af is 


NRfe 




k 


rV 


k 


Pi “[«* +n/vi'2 


rye 


( 10 ) 


where N(M) is the number of nearest (next)-neighbors in 
the equivalent crystal of species i and R { (R 2 ) is the 
nearest (next)-neighbor distance in the equivalent crystal 
of lattice parameter af. r l and r 2 , are the equilibrium 
nearest- and next-nearest-neighbor distances in an etui- 
librium crystal of species /, respectively. The chemical 
energy is then computed with 

Ef=yE , c F*{af*) 


( 11 ) 
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TABLE I. Experimental input: Cohesive energy (in eV), lattice parameter (in A). Equivalent crystal 
theory parameters (Ref. 11 ): p,l (in A), a (in A '), and k (in A) for several fee elements. The surface en- 
ergy (in ergs/cm 2 ), computed with ECT (Ref. 1 1), is also included. The last entry displays the BFS pa- 
rameters and A fl/I for the Ni-Au system. 


Element 

Cohesive 

energy 

Lattice 

constant 

P 

/ 

a 

X 

Surface 

energy 

Au 

3.78 

4.078 

10 

0.236 

4.339 

0.663 

1621.55 

Ni 

4.435 

3.524 

6 

0.270 

3.015 

0.759 

3076.26 



BFS: A Au Ni “ 

-0.050 62, A NiAu = -0.062 25 




and 

$=YoEc F *(af 0 *) , ( 12 ) 

where y(y 0 ) — + 1 if 0/ C *(fl, ° ) - 0 and Y^Yo^ = 1 °th _ 
erwise, and af* =q(af~a ^ )//,. The scaled lattice pa- 
rameter af* is obtained from Eq. (10) with the parame- 
ters a ik listed in Ref. 8, and a% * is computed by solving 
Eq. (10) but with a ik = a,. The rest of the parameters ap- 
pearing in Eq. (10) are listed in Ref. 1 1. 

Even though BFS is a semiempirical method, its depen- 
dence on experimental input is minimal in that only two 
experimental (or theoretical) alloy values (in the present 
study the heats of solution in the dilute limit were used 12 ) 
are needed. The remaining input are pure element prop- 
erties: the cohesive energy, equilibrium bulk modulus, 
and lattice parameters. In this work, we used the param- 
eters A^ and & BA determined following the procedure 
outlined in Ref. 8. The experimental input, as well as the 
resulting BFS parameters can also be found in Ref. 8. 
The BFS and ECT parameters used in this work for Ni 
and Au are listed in Table I. 

Before proceeding with the application of BFS to the 
problem of surface alloying, we should emphasize that in 
the context of BFS, the strain and chemical energy con- 
tributions differ substantially in meaning from the one 
these terms have in other approaches. The BFS strain 
energy is related to the usual strain only in that the atom- 
ic locations are those found in the actual alloy: the BFS 
strain energy of a given atom is then the actual strain 
that it would have in a monatomic crystal of the same 
species of the reference atom. Likewise, the BFS chemi- 
cal contribution is related to the usual chemical energy in 
that the actual chemical composition of the alloy is taken 
into account, but with the neighboring atoms located in 
ideal atomic sites: the BFS chemical energy of a given 
atom is then the actual chemical energy in an ordered en- 
vironment with the lattice spacing characteristic of the 
equilibrium lattice of the reference atom. Therefore, the 
BFS contributions are, in a sense, a certain combination 
of the actual strain and chemical energies. We refer the 
reader to previous applications of BFS for more insight in 
this issue. 8,9 

III. RESULTS AND DISCUSSION 

In this BFS study of surface alloying, we will concen- 
trate on the low-coverage regime, where the alloying 


effect is more pronounced. 1 Two approximations were 
made: while the STM experiments were carried out at 
room temperature, 1 our simulation was done at zero tem- 
perature, in the belief that the essential ingredients that 
drive the alloying process will be present at zero tempera- 
ture. Also, for the sake of simplicity, we ignored lattice 
relaxations due to either the surface or the presence of 
adatoms with large lattice mismatch with the substrate. 
This simplification allows for the study of large numbers 
of possible configurations without much computational 
effort. Moreover, previous effective-medium theory cal- 
culations, 1 also done at zero temperature, indicate that 
both approximations are reasonable in the sense that evi- 
dence for surface alloying is found even if relaxation is 
left out of the simulation. 

The EMT calculation 1 succeeds in explaining the basic 
feature of this phenomenon, namely, the energy gain real- 
ized when Au atoms substitute Ni atoms only in the sur- 
face plane. Arguments based on the cohesive energy 
function of each element clearly substantiate this numeri- 
cal finding. However, the EMT calculation predicts the 
exact same energy gain for the substitution of Au mono- 
mers and dimers, thus failing to explain the abundance of 
dimers in the experimental situation. Moreover, some 
disagreement exists with regard to the location of the 
ejected Ni atoms. While it would be highly unlikely that 
a theoretical model — such as the EMT calculation or the 
one in the present work — could reproduce ail the details 
of the experiment, it is a reasonable demand to expect to 
find evidence for the most salient features among the pre- 
dicted outcomes. With that in mind, we will direct this 
application of BFS to answer the question if it is possible 
to obtain the types of configurations that are inferred 
from experiment. 

There are at least two ways of performing this simula- 
tion: a systematic one, where Monte Carlo techniques can 
be used to determine the equilibrium configurations for a 
certain coverage, or a brute-force one, where specific 
configurations are chosen to study the detailed processes 
taking place and the behavior of individual atoms. We 
chose the latter: by examining a sufficiently large number 
of possible configurations, even those that are energetical- 
ly unfavorable, we expect to gain some insight that could 
later be applied for similar systems and therefore search 
for a general criterion. 4 

The calculation was performed on a Ni slab several 
layers deep, with a (1 10) surface. Varying numbers of Au 
atoms were deposited and located in substitutional sites 
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in the top or inner layers, or just as adatoms on hollow 
sites of the Ni(l 10) surface. 

We will analyze the results in terms of the energy of 
formation of a given configuration, and the contributions 
of individual atoms, as defined in Eq. (3), to that magni- 
tude. Let 8 H be the energy of formation per impurity 
atom (in eV/atom): 


bH — 


Atf-A H 0 

*Au 


(13) 


where A H is the energy of formation of the configuration, 
A H q is the corresponding value for a free Ni surface and 
N Au is the number of impurity atoms. Starting with a 
single Au adatom deposited on a hollow site on the Ni 
substrate, it is seen that it contributes so that the energy 
of formation of the cell considered is reduced by —0.629 
eV. Within the framework of BFS, this is due to the de- 
crease of the BFS strain energy of the four Ni atoms on 
the surface plane as they increase their coordination by 
having the Au atom as a nearest neighbor. Furthermore, 
the BFS chemical energy contributions from these Ni 
atoms and the Au atom are also negative. The change in 
energy due to the addition of a second Au adatom de- 
pends on its location in the overlayer with respect to the 
first: 6H is —0.629 eV if the two adatoms are far from 
each other, and it varies from —0.624 to —1.08 eV for 
dimers located in the diagonal direction and the close- 
packed (cp) direction, respectively. Additional adatoms 
have basically the same effect: clustering along the cp 
direction always reduces bH , with the formation of a Au 
chain being the preferred configuration. As we will see in 
detail later for JV Au =4, it is interesting to note that two 
configurations involving Au dimers are very close in ener- 
gy to the Au chain, always along the cp direction. This 
supports the claim of the likelihood of dimers being 
formed on the overlayer before the exchange with Ni 
atoms begins. 

Returning to the case of just one Au atom, it is in- 
teresting to study the different configurations in terms of 
the strain, glue, and chemical energy, as described by Eq. 
(3): bH results from a delicate balance between these 
quantities. It was noted above that while the chemical 
energy contribution from the Ni-Au bonds is negative, its 
effect in lowering bH is modulated by the value of the 
glue term, g, which in turn depends on the magnitude of 
the strain energy contribution e 5 [see Ref. 8 for a detailed 
description of the calculation of the different terms in Eq. 
(3)]. As the impurity atom penetrates into the Ni sub- 
strate occupying sites on the surface plane and the planes 
below, it reduces its BFS strain energy significantly, as it 
finds increasing coordination as well as a much higher 
electron density due to the difference in size between Au 
and Ni atoms. Therefore, with e Au lower, the glue term 
increases in magnitude, emphasizing the negative contri- 
bution to bH due to the chemical energy. However, 
beyond the surface plane the BFS strain energy increases 
again as the Au atom finds itself in a compressed bulklike 
environment, thus reducing the glue and the negative 
contribution of the chemical energy. A quick estimate il- 
lustrates this argument: the calculation of the BFS strain 


energy (to a nearest-neighbor approximation) is basec on 
a measure of the defect as seen by a given atom. Equa- 
tion (9) establishes a relationship between the defect crys- 
tal (right-hand side of the equation) and the equivalent 
crystal (left-hand side). The term to the left could be un- 
derstood as a measure of the defect, given by (ii a 
nearest-neighbor approximation) q d =nr p e ~ ar , 8,9 where n 
is the number of nearest neighbors located at a distan:e r 
of the atom in question (assuming, as is the case in his 
unrelaxed calculation, that the separation distance be- 
tween nearest neighbors is the same in all cases). The pa- 
rameters p and a depend on the species of the reference 
atom. 8,9,11 Equilibrium [a situation for which both sides 
of Eq. (8) are identical^ corresponds to q e ~Nr p e ar % 
where N = 12 and r e —V 2 a e /2 (for fee metals). For a Au 
atom in a Ni lattice, we could ask ourselves what the 
effective number of Ni nearest neighbors that will simu- 
late the equilibrium situation for a Au atom (i.e., what is 
the value of n for which q d =q e ) is 

nr^ u e aAu ' Ni = iVr A * u e y < 14) 

where r x is the equilibrium nearest-neighbor distance for 
species X , p Au ^10, and a Au = 4.339. n It turns out that 
n ~9.4, which means that if the Au atom is surrounded 
by n Ni atoms in a Ni lattice it would have no BFS stiain 
energy. Conversely, a Au atom in a substitutional site in 
the Ni(110) surface with seven Ni nearest neighbors at Ni 
equilibrium distances has the same strain energy that it 
would have in a Au lattice with 8.9 Au nearest neighbors. 
If the number of nearest neighbors was a continuous vari- 
able, the Au atom would be in equilibrium somewhere be- 
tween the top two layers of the Ni slab. 

Table II displays the values of bH for the Au atom lo- 
cated (a) in the overlayer (O), (b) in a substituted site in 
the surface plane (S) with the substituted Ni atom in the 
overlayer as a nearest neighbor, (c) same, with the substi- 
tuted Ni atom in the overlayer far from the impurity, (d) 
in the first plane below the surface (16), and (e) two 
planes below the surface plane (2b). The intermediate 
columns indicate the values of e Au , g Au , and e Au — : A ° 

( bH , of course, includes the contributions from the sur- 
rounding Ni atoms). 

The small difference in bH between configurations (b) 


TABLE II. One Au atom in different locations (see text). 
The strain energy, glue, and chemical energy contributions are 
listed in the second, third, and fourth column, respectively. The 
fifth column displays the total contribution to the energy of or- 
mation from the Au atom and the last column shows the total 
energy of formation per impurity atom of the cell. All energies 
are in eV. 


Config. 

*^Au 

8 Au 

_ c o 

b Au fc Au 

^Au 

bH 

a: ( O ) 

1.3670 

0.2810 

-0.3506 

1.2683 

-0.62V 12 

b: (S) 

0.2125 

0.6840 

-0.4574 

-0.1003 

-0.761 45 

c: (S) 

0.5326 

0.5186 

-0.4469 

0.3009 

-0.79: 18 

d: (16) 

0.6873 

1.6610 

-0.2449 

0.2806 

0.665 55 

e: (26) 

3.5988 

2.6699 

-0.9820 

0.9581 

1.475 55 
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TABLE III. Contributions per layer to the energy of formation per atom (in eV) of a free Ni surface 
(a), a surface with a Au atom in a substitutional site in the surface plane, with the substituted Ni atom 
nearby (b), and the same case when the Ni atom is somewhere else in the overlayer (c). 


Configuration 
(see text) 

A// 

overlayer 

A H 

surface 

AH 

one-below 

AH 

two-below 

b H 

(a) Free surface 

0.0 

1.6097 

0.13198 

0.0 


(b) 

2.4516 

1.4093 

0.11835 

-0.0002 

-0.762 45 

(c) 

2.4604 

1.4080 

0.11703 

-0.0002 

-0.792 18 


and (c) can be easily explained in terms of the larger num- 
ber of surface Ni atoms in (c) affected by the impurity 
atom and the substituted Ni atom in the overlayer: this 
is illustrated in Table III, where the contribution from 
each layer is listed for each case. As a reference, we also 
indicate the corresponding contributions from a free Ni 
surface (i.e., no Au atoms present). 

Most configurations with two impurity atoms are ener- 
getically favored with respect to those with a single Au 
atom: our results show that a Au dimer immersed in the 
surface layer, with the substituted Ni atoms forming a di- 
mer somewhere else in the overlayer (i.e., not sharing 
next-nearest neighbors) have the lowest energy (with both 
dimers oriented in the cp direction) with respect to the 
two Au atoms and two substituted Ni atoms in other lo- 
cations: following the structure of Table III, Table IV 
displays the contributions per atom from different layers 
for the following configurations: (a) two isolated Au ada- 
toms in the overlayer [(0) two adatoms], (b) a gold dimer 
in the overlayer [(0) Au 2 ], (c) a gold dimer in the surface 
plane with a Ni dimer in the overlayer, nearest neighbors 
[(S)Au 2 ,(0)Ni 2 ,NN], (d) a gold dimer in the surface 
plane with a Ni dimer somewhere else in the overlayer 
[(S)Au 2 ,(0)Ni 2 , far], (e) a gold dimer one plane below the 
surface with the Ni dimer somewhere else in the over- 
layer [( 16)Au 2 ,(0)Ni 2 , far], and (f) a gold dimer two 
planes below the surface with the substituted Ni dimer in 
the overlayer [(26)Au 2 ,(0)Ni 2 ]. These last two cases are 
included to highlight the fact that Au dimers penetrate, 
at the most, into the surface layer. We also list the ener- 
gy of formation of the cell per impurity atom as well as 
the results for a free surface. The last entry in Table II 
and the last entry in Table IV show that, in spite of the 


ordering found in the surface plane (i.e., the tendency of 
Au atoms to form dimers and occupy substitutional sites 
in that plane), the Ni-Au system phase separates in the 
bulk. 

These results underscore the possibility that after depo- 
sition dimers tend to form on the overlayer along the 
close-packed direction, and later occupy substitutional 
sites in the Ni substrate, but only in the surface plane, a 
characteristic feature observed in the STM images. The 
tendency for dimer formation in the adlayer is also sup- 
ported by the EMT results. 

The results in Tables II and IV indicate that two isolat- 
ed Au adatoms realize a greater gain in energy by 
diffusing along the surface and forming a dimer than if 
they exchange places with Ni surface atoms. With the 
diffusing mechanism becoming more important at higher 
temperatures, one would then expect dimer formation in 
the adlayer prior to exchange to be favored. The final 
state, [(S)Au 2 ,(0)Ni 2 , far], that can be reached if the Au 
diffusion mechanism dominates, is however higher in en- 
ergy than another alternative, [ j (S)Au) 2 ,(0)Ni 2 , far], a 
state that can only be reached if the exchange mechanism 
of isolated adatoms dominates. The EMT results, while 
indicating a strong energy gain in dimer formation, favor 
Au insertion in the surface over Au ad-dimer formation 
and yield no difference between the adsorbed monomer or 
dimers, which, as noted in Ref. 1, fails to explain the 
abundance of dimers observed experimentally. The re- 
verse is true in the BFS scheme: Au dimer formation is 
strongly favored with respect to isolated Au substitutions 
therefore increasing the likelihood of finding the dimer 
substitution process over monomers. Whether this argu- 
ment correctly describes the actual mechanism should be 


TABLE IV. Contribution per layer to the energy of formation per atom bH indicates the energy of 
formation per impurity atom for several configurations with two Au atoms (0.036 ML). 

Configuration 
(see text) 

AH 

overlayer 

AH 

surface 

AH 

one-below 

AH 

two-below 

bH 

Free surface 

0.0 

1.6097 

0.131 98 

0.0 


(O) Two adatoms 

1.268 30 

1.4141 

0.11679 

0.0 

-0.629 12 

(0) Au 2 

0.806 82 

1.4151 

0.11687 

0.0 

-1.080 63 

(S)Au 2 ,(£»Ni 2 ,NN 

2.019 53 

1.2915 

0.11016 

-0.0004 

- 1 .043 67 

(S)Au 2 ,( 0)Ni 2 , far 

2.037 69 

1.2812 

0.10741 

-0.0004 

- 1.14302 

((S)Au),,(0)Ni 2 , far 

2.037 69 

1.5096 

0.124 53 

-0.0001 

- 1.192 94 

(lf>)Au 2 ,(0)Ni 2 

2.037 69 

1.4176 

0.14167 

-0.0022 

0.373 05 

(2b)Au 2 ,(0)Ni 2 

2.037 69 

1.4229 

0.108 52 

-0.1701 

1.653 13 
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further corroborated by detailed first-principles calcula- 
tions. To the effect of this work, where we are trying to 
simulate the experimental evidence with a theoretical 
model, this analysis is just meant to provide a measure 
for the interpretation of results corresponding to higher 
coverages. 

Obviously, increasing the number of deposited Au 
atoms leads to countless configurations, impossible to an- 
alyze in any systematic fashion. With our goal being to 
determine the likelihood of the experimental interpreta- 
tion to have a theoretical match, we, therefore, limit the 
following examples to those cases where every multiatom 
arrangement is restricted to be oriented in the cp direc- 
tion and the gold atoms are restricted to be located only 
in the overlayer and/or the surface plane. In a cell with 
60 atoms in each plane, the case with 7V All = 4 corre- 
sponds to a coverage of 0.067 ML. Assuming periodicity, 
the cp chain [Fig. 1(a)] has the lowest energy per impuri- 
ty atom among all the possible configurations with all 
four Au atoms in the overlayer (5// = — 1.290 eV) as 
compared to other configurations shown. In Fig. 1, Ni 
atoms are indicated by small circles (surface plane) or 
large circles (overlayer), Au atoms are correspondingly 
indicated by small and large disks. Also, for reasons of 
space, we limit all the figures in this work to display only 
the active region of the 60 atoms cell (i.e., the region 
affected by the Au atoms and the substituted Ni atoms). 
Not surprisingly, the next possible configuration corre- 
sponds to two cp dimers, far from each other [Fig. 1(b)]. 
This is followed by different island shapes, as illustrated 
in Figs. 1(c)- 1(h). The corresponding values of 6 H are 



FIG. 1. Schematic representation of the Ni(llO) surface: 
small circles (o) indicate Ni atoms in the surface layer; Au 
atoms in the overlayer are indicated by large disks (O). These 
configurations correspond to a Au coverage of 0.067 ML. 


TABLE V. Energy of formation per impurity atom 8H for 
Au coverage of 0.067 ML for the configurations indicated in 
Fig. 1. 


Config. 
(Fig. 1) 

8H 

(eV/atom) 

Config. 
(Fig. 1) 

6 h 

(eV/atom) 

(a) 

-1.29000 

(e) 

-1.062 48 

(b) 

-1.08063 

(f) 

-1.061 91 

(c) 

-1.078 30 

(g) 

-1.053 11 

(d) 

-1.07193 

(h) 

-1.037 37 


listed in Table V in order of decreasing energy. 

However, there are several configurations with louver 
energy for the same coverage: these correspond to the 
case when the Au atoms substitute for Ni atoms in the 
surface layer, with the displaced Ni atoms forming cp 
chains of four atoms in the overlayer. The difference be- 
tween these configurations is in the relative positior of 
the Au atoms inserted in the surface layer. Figure 2 indi- 
cates six of the lowest-energy configurations and Table VI 
lists the corresponding energies: all of them are lower 
than (a) in Table V (the lowest-energy configuration w ith 
all the Au atoms in the overlayer), indicating that ar- 
rangements that include long Ni chains surrounded by 
Au dimers are among the most likely configurations. " r he 
lowest-energy configuration [Fig. 2(a)] is characterized by 
the insertion of individual Au atoms and the formation of 
a Ni chain, while all the other configurations [except Fig. 
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FIG. 2. Schematic representation of the Ni(l 10) surface: as in 
Fig. 1, small circles (o ) indicate Ni atoms in the surface layer; 
large circles (O) indicate Ni atoms displaced to the overlayer; 
Au atoms are indicated by large disks (O) when in the overlayer 
and with small disks (o ) when occupying Ni sites in the surface 
plane. These configurations also correspond to a Au coverage 
of 0.067 ML. 
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TABLE VI. Energy of formation per impurity atom S H for 
Au coverage of 0.067 ML for the configurations indicated in 
Fig. 2. 


Config. 
(Fig. 2) 

bH 

(eV/atom) 

Config. 
(Fig. 2) 

bH 

(eV/atom) 

(a) 

-1.39647 

(d) 

-1.342 86 

(b) 

-1.348 72 

(e) 

-1.341 22 

(c) 

-1.346 64 

(0 

- 1.335 80 


2(d)] display Au dimers in different locations relative to 
the Ni chain. The spread in energy between these states 
is quite small (0.06 eV). This fact, together with the ap- 
proximations made in this calculation, raises a question 
with regard to which one is the true ground state. In 
spite of this, an argument could be made with respect to 
the actual process. The fact that Fig. 2(a) has the lowest 
energy suggests that the insertion of isolated Au atoms is 
preferred, but as discussed above, the mechanisms lead- 
ing to this final state might be less favored than those 
leading to alternative configurations, i.e., the highly, en- 
ergetically favored formation of adlayer dimers over in- 
sertion at elevated temperatures. 

From these results one can see, even at this very low 
coverage, some indication of the trends which ultimately 
would lead to the situations found experimentally. These 
configurations share some distinctive features: the 

penetration of Au atoms in the surface layer, the forma- 
tion of Ni chains in the overlayer along the close-packed 
direction, and the linkage of the Au atoms and the substi- 
tuted Ni atoms by means of an intermediate surface Ni 
atom. We illustrate these last two issues in Fig. 3. The 
distinction between Figs. 3(a) and 3(b) is given by the 
shape of the Ni island: Fig. 3(a) shows a four-atom Ni 
chain while Fig. 3(b) displays the same chain but broken 
into two separate Ni dimers, with a substantial increase 
in energy as indicated in Table VII. 

An interesting detail is the fact that the configuration 
shown in Fig. 3(a) has a lower energy bH than the one 
shown in Fig. 3(f), the difference between the two being 
the orientation of the Au dimers in the surface layer. 
Our previous discussion would lead us to expect the 
second one to be lower in energy because of the orienta- 
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FIG. 3. Using the convention indicated in Fig. 2, the 
configurations shown correspond to a Au coverage of 0.067 ML. 


tion of the dimer in the close-packed direction, but that is 
not the case. On the one hand, it is clear that the advan- 
tage for the Fig. 3(a) configuration is that it maximizes 
the gain in energy due to the particular linkage between 
the Au atoms and the substituted Ni atoms: by position- 
ing themselves perpendicular to the close-packed direc- 
tion, all four Au atoms benefit from the gain due to the 
linkage effect discussed above. On the other hand, the 
question arises if the configuration shown in Fig. 3(a) has 
a high probability for existence. The answer is clearly no: 
as argued earlier, there is a strong likelihood that the Au 
atoms deposited on the Ni substrate will migrate to form 
dimers oriented in the close-packed direction. While 
6H= — 1.0806 eV for a Au dimer in the close-packed 


TABLE VII. Contribution per layer to the energy of formation per atom. bH indicates the energy of 
formation per impurity atom for several configurations with (a) and (b) four Au atoms (0.067 ML), 
(c)-(e) a Au and a Ni dimer (0.033 ML), and (f)-(h) Au dimers and a Ni chain (0.067 ML). 


Configuration 
(see text) 

AH 

overlayer 

A H 

surface 

A H 

one-below 

AH 

two-below 

bH 

Fig. 3(a) 

1.8251 

1.4097 

0.1173 

-0.0002 

- 1.398 25 

Fig. 3(b) 

2.0377 

1.4092 

0.1172 

-0.0002 

-1.194 64 

Fig. 3(c) 

1.0188 

1.5111 

0.1246 

-0.0001 

- 1.14302 

Fig. 3(d) 

1.0188 

1.5110 

0.1246 

-0.0001 

-1.145 12 

Fig. 3(e) 

1.0158 

1.5114 

0.1249 

-0.0001 

-1.132 18 

Fig. 3(0 

1.8251 

1.4132 

0.1173 

-0.0002 

-1.346 64 

Fig. 3(g) 

1.8251 

1.4130 

0.1173 

-0.0002 

-1.348 72 

Fig. 3(h) 

1.8191 

1.4136 

0.1179 

-0.0002 

- 1.335 80 
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direction, it is just —0.6287 eV for a dimer in the perpen- 
dicular direction. Moreover, two isolated Au adatoms 
have an even lower energy than that (bH = —0.6291 eV). 
This would be followed by an exchange of Au and Ni di- 
mers, with the substituted Ni dimers forming islands in 
the overlayer. It was shown that, once again, dimers will 
conserve their orientation after the exchange process. Fi- 
nally, once squeezed out from the surface layer, the Ni 
atoms will migrate on the substrate forming islands with 
close-packed chains as their basic structure. Therefore, 
in an experimental case, the configurations shown in Figs. 
3(a) and 3(b) are highly unlikely to be found. 

The remainder of Fig. 3 shows three rather similar 
configurations for two different Au coverages (0.033 and 
0.067 ML), where the difference resides in the relative lo- 
cation of the Au dimers and the Ni chain: the gain in en- 
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ergy is maximized when they are connected by a Ni sub- 
strate atom [boxed in Figs. 3(d) and 3(g)]. Table VII lists 
the energy contribution per layer, showing the small 
changes taking place when evolving from Figs. 3(c) to 3(e) 
(N Au =2) and Figs. 3(f) to 3(h) (tf Au =4). 

As we increase the coverage, a pattern emerges: the en- 
ergy spectrum of the large number of configurations 
available shows a tendency to group into energy bands in 
the sense that within each group there are small 
differences in energy, while the energy gap between each 
band is much larger. Also, each group of configurations 
are characterized by a certain symmetry: in the lower end 
of the spectrum, we always find configurations where the 
essential features are Au atoms in the surface plane, 
forming dimers, and the substituted Ni atoms for rung a 
chain along the cp direction in the overlayer, as shown in 
Fig. 4. The difference among the configurations belong- 
ing to this group is given by the location of the dimers 
with respect to the Ni chain. These findings aie con- 
sistent with the large dimer concentration seen in the 
STM experiments. 1 

However, as the coverage increases, configurations ap- 
pear in this ground-state group: those where the Au ada- 
toms form long chains along the cp direction in th? over- 
layer [Fig. 4(d)]. This could be taken as an indication of 
growth of the Au film on the substrate in competition 
with the formation of a surface alloy. A first hint of this 
alternative can already be seen at relatively low cover- 
ages: Table VIII shows some results for 0.13- and 0.17- 
ML Au coverage, with the corresponding configu ations 
represented in Fig. 4. 

A possible explanation for this change in grow h pat- 
tern can be found in terms of the surface energy of Au be- 
ing much lower than that of Ni (see Table I). For low 
coverages, the decrease in energy driving the penetration 
of Au atoms in the surface layer is guided by the effective 
coordination effect mentioned earlier: Au atoms benefit 
from locating themselves in the surface layer, with the Ni 
atoms forming islands in the overlayer. At one point, the 
increase in surface energy due to the large Ni islands be- 
comes larger than any gain generated by the inter nixing 
of Au and Ni atoms in the surface plane, therefore, 
configurations with Au islands on the Ni substrate be- 
come energetically favored thus reverting to a normal 
growth mode where Au atoms tend to form a pure Au 
layer. The breaking point between these two regimes 
seems to be around a Au coverage of 0.5 ML. For higher 
coverages, there is experimental indication that alterna- 
tive 3D pattern formation starts. This will be the subject 
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FIG. 4. Using the convention indicated in Fig. 2, the 
configurations shown correspond to a Au coverage of 0.13 ML 
[for (a), (b), and (c)] and 0.17 ML [for (d), (e), and (f)]. 


TABLE VIII. Energy of formation per impurity atom bH for 
Au coverage of 0.13 ML (left columns) and 0.17 ML (right 
columns) f or the configurations indicated in Fig. 4. 


Config. 
(Fig. 4) 

bH 

(eV/atom) 

Config. 
(Fig. 4) 

bH 

(eV/atom) 

(a) 

-1.44847 

(d) 

- 1.497 35 

(b) 

- 1.393 17 

(e) 

- 1.393 89 

(c) 

- 1.335 80 

(f) 

- 1.378 70 
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TABLE IX. Contributions per layer to the energy of formation per atom (in eV) of a free Ni surface, 
and each of the three configurations shown in Fig. 5 (Au coverage of 0.50 ML). 

Configuration 
(see text) 

Ai/ 

overlayer 

A H 

surface 

A H 

one-below 

A H 

two-below 

8 H 

Free surface 

0.0 

1.6097 

0.13198 

0.0 


Fig. 5(a) 

1.8228 

0.2284 

0.033 97 

-0.001 63 

-1.13900 

Fig. 5(b) 

1.5134 

0.8090 

0.060 50 

0.000 75 

-0.229 43 

Fig. 5(c) 

0.4823 

0.7866 

0.083 53 

0.002 17 

-1.25644 


of future work. To illustrate these observations, we con- 
sider a few configurations corresponding to 0.5-ML cov- 
erage. Figure 5 shows three configurations with 30-Au 
atoms in a 60 atoms-per-plane cell. Figure 5(a) shows a 
highly disordered distribution, dominated by the pres- 
ence of Au dimers inserted in the surface plane, separated 
by irregular Ni islands in the overlayer. All 30 Au atoms 
are located in the surface. Figure 5(b) shows a highly or- 
dered distribution, where the Au atoms are sandwiched 
between two Ni layers. Figure 5(c) displays a large Au is- 
land on the pure Ni substrate. The corresponding nu- 
merical results are listed in Table IX: the contribution of 
each layer (following the format in Table II) to the total 
energy of formation is shown, confirming our previous as- 
sumption. For completeness, we also show results for the 
pure Ni slab (no Au coverage) to highlight the surface 
effects generated by the presence of the Au atoms and 
their distribution. The dimer + island configuration [Fig. 
5(a)] is characterized by Au dimers loosely linked to the 
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FIG. 5. Using the convention indicated in Fig. 2, the 
configurations shown correspond to a Au coverage of 0.50 ML. 


Ni islands: if the Au dimers where covered by substituted 
Ni atoms, as is the case in Fig. 5(b), they would have their 
strain energy increased to levels where little gain is real- 
ized from the substitution process: the contribution from 
the surface layer jumps from 0.2284 eV in Fig. 5(a) to 
0.8090 in Fig. 5(b). Also, in Fig. 5(b) there is a large con- 
tribution from the Ni overlayer (1.5134 eV) due to the 
high surface energy of Ni. In Fig. 5(c), the reversal 
brings stabilization to a sandwich distribution: the Au 
overlayer has a noticeable lower contribution (0.4823 eV 
vs 1.5134 for Ni) and the surface of Ni lowers its energy 
substantially (0.7866 eV) with respect to a free Ni surface 
(1.6097 eV). 


IV. CONCLUSIONS 

Summarizing the results presented for different cover- 
ages, we can imagine the sequence of events that could 
lead to the situation observed experimentally: The first 
few Au adatoms are readily adsorbed on the Ni surface 
with a tendency at elevated temperatures to migrate to 
form the energetically favored dimers along the close- 
packed direction. These dimers are then exchanged with 
Ni atoms on the surface plane conserving their original 
orientation. The displaced Ni atoms tend to form islands 
along the close-packed direction trying to keep a certain 
level of linkage with the Au dimers embedded in the sur- 
face layer. This process leads to the formation of islands, 
whose ultimate shape is, therefore, determined by the rel- 
ative location of the Au dimers in the surface plane. The 
low energy of formation of the configuration shown in 
Fig. 2(a) suggests that isolated Au atoms inserted in the 
surface plane are also likely to be found. For low cover- 
ages (less than 0.5 ML), these arguments provide a plausi- 
ble explanation for the observed experimental results. 
Two processes compete to bring a delicate balance that 
essentially favors the formation of a surface alloy: the en- 
ergetically favorable intermixing of Au atoms in the sur- 
face plane due to the increased effective coordination per- 
ceived by those surface atoms, and the energetically un- 
favorable formation of islands with the substituted Ni 
atoms, with an increase in energy due to the lower coor- 
dination. As the size of these islands grows, more energy 
gain is realized by a direct deposition of a Au overlayer 
as opposed to the formation of a surface alloy. The final 
configuration is, obviously, strongly dependent on the ex- 
perimental conditions, as it could also be possible that in 
a slow deposition process, the Ni islands could allow for 
the formation of additional surface alloy cells, therefore 
generating a disordered alloy pattern that goes beyond 



7216 


BOZZOLO, IBANEZ-MEIER, AND FERRANTE 


51 


the single-layer alloy formed at low coverages. A sudden 
coating of Au would not allow for the formation of any 
surface alloy at all. 

The experimental data available for high coverages 
shows evidence of interesting growth patterns, which we 
plane to analyze in future efforts: undoubtedly, these 
features are a direct consequence of the low miscibility 
between the two participating elements and a careful nu- 


merical study of the type presented in this work might be 
helpful to gain insight into this growth process. 
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